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ABSTRACT 


This thesis analyzes the performance of a ramjet with an 
Oswatisch inlet using a blunt centerbody and compares perfor- 
mance to a baseline ramjet using an inlet with a conical 
Spike at Mach 3.0. 

Inlet performance as a ratio of inlet lip to nose center- 
body ratio, r,/T,> is developed. The capture streamline for 
each ratio is determined and the coefficient of additive drag 
is calculated as a function of r,/Th: Setting thrust coef- 
ficient equal to coefficient of drag, the performance of two 
ramjets is determined. One ramjet is the baseline with a 
Spake inlet; the other ramjet uses the blunt centerbody. 

Ramjets and inlets are compared on the basis of specific 
fuel consumption, excess thrust coefficient and specific 
thrust. For the ramjet with blunt centerbody, performance 
parameters were calculated as a function of inlet lip radius 
to nose centerbody radius. Also compared is the effect of 
the ratio, r,/Tp> on relative detection range. For both 
types of ramjets, the detection range is reduced by approxi- 
mately 66%. Performance of the ramjet with blunt nosed 
centerbody is severly handicapped due to high additive 
drag and poor pressure recovery. Specific fuel consumption 
1s approximately 50% greater for the ramjet with the blunt 


centerbody compared to the ramjet with the spike inlet. 
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io LNERODUCTION 


A. BACKGROUND 

The rapid development and operational deployment of long 
range antiship cruise missiles has generated demanding 
requirements for the anti-ship missile defense (ASMD) mis- 
Sion. The threat has necessitated that the Navy develop the 
concept of a "defense in depth" for the defensive doctrine 
of a carrier battle group (CBG). The doctrine emphasizes 
the utilization of all available weapon systems in a layered 
defense to defeat hostile targets. Conceivably, an anti- 
Ship missile (ASM) could survive attacks from extremely long 
range missiles, the E-2C and F-14 combat air patrol (CAP) 
team, Standard ER and Standard MR missiles, conventional gun 
ordnance, point defense weapons and close in defense weapons, 
such as Phalanx, prior to reaching a high value unit. Since 
enemy doctrine is to fire the ASM in such numbers as to 
overwhelm the Command and Control capability of the CBG 
and of an individual vessel, the accent of new ASMD weapons 
Systems is on quick response time, short time of flight and 
a high degree of accuracy against a manuvering ASM. 

Although a single saturation raid could possibly be 
defeated, what about a second or third raid in a single day? 


Will the CBG have enough missiles to supplement airborne CAP 
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and still destroy the ASM threat 30 to 40 miles from the 
et, Value unit? 

One solution to this vexing problem is to modify the 
5"/54 Semi-Active Laser Guided Projectile (SALGP) which is 
[melistrated if Figure I[-l1. The SALGP round is currently un- 
dergoing operational evaluation (OPEVAL) prior to fleet wide 
gmtroduction [Ref. 1] and consists of a five inch projectile 
with a laser seeker on the projectile nose. The round is 
designed to “home in" on laser energy reflected from a target 
by a laser illuminator. The modification to the SALGP round 
consists of replacing the laser seeker head with an infrared 
(IR) seeker and adding a solid fuel ramjet to power the pro- 
jyectile. The ramjet would provide constant speed over a 
longer portion of the flight path and provide energy in the 
terminal phase of flight to allow for manuvering to hit a 
manuvering target. The modified SALGP is identified as a 
gun launched missile (GLM). The modifications would allow 
the GLM to become a fire and forget round requiring no ter- 
minal guidance from the ship. 

Although not helping in the long range defense of a CBG, 
the GLM would complement the point defense and close in 
defense weapon systems to the ships in the U.S. Navy today. 

The GLM will allow a naval vessel without a missile sys- 
tem and with a properly configured gun system to increase 


its anti-air warfare (AAW) capability in a cost effective 
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way. More GLM rounds can be carried in existing gun maga- 
zines than missile rounds can be carried in current missile 
magazines. Within a ten nautical mile zone around a ship, 
the GLM will have a faster response time, a shorter time af 
flight and be as accurate as guided missiles in the current 
Navy inventory. Fast response and short time of flight are 
due to the initial velocity imparted to the GLM by the 5'/54 
gun. The shape of the nose of the GLM round will be con- 
strained by the use of the hemispherical lense surface for 
the infrared seeker in the projectile nose. This blunt nose, 
although optimized for the sensor optics, is not necessarily 
the best aerodynamic shape. 

Although the SALGP has the potential to be highly effec- 
tive in the Naval Gunfire Support (NGFS) and in the Surface 
to Surface Warfare (SUW) missions, it was not designed ini- 
tially for the Air Warfare (AAW) or Anti-ship Missile 
Defense (ASMD) missions. Consider a GLM based on SALGP 
hardware to the maximum extent possible. To become success- 
ful in the AAW/ASMD mission areas, it is essential to im- 
prove the overall SALGP performance and capability. One area 
to achieve easily measured performance gains is in the 
method of propelling the SALGP. 

One propulsion proposal that decreases mission flight 
time and increases mission effectiveness is to replace the 


TMG Eemeillant rocket on the SALGP round with that of a 


18 





pameeewttieettier Jiquid or solid fuel. Obviously, to keep 
total system developement and procurement costs low, and 
therefore warrant the use of a ramjet in a cost effective 
analysis, the SALGP round should be changed as little as 
Gomcable in adding the ramjet. Effectively, the designer is 
constrained to use a blunt nosed centerbody to house the IR 
seeker and sensor package. The remainder of the projectile 
is constrained by the dimensions of the handling equipment 
in the MK 45 mount and Department of the Navy Specifications 


[Ref. 2]. 


Be PURPOSE 

One purpose of this thesis is to investigate the effect 
that a blunt nosed centerbody has on the pressure recovery 
Seereintet. lhe pressure recovery iS an important indica- 
tor of ramjet performance, for poor pressure recovery will 
reduce the flight envelope over which the ramjet operates. 

Although drag is large, the blunt nosed body was aerody- 
Namically feasible for a rocket guided projectile, since the 
operation of the rocket engine does not require air to be 
"brought onboard". The use of a ramjet may prove inadvise- 
able or infeasible due to aerodynamics of the flow around 
the blunt nose body enroute to the inlet annulus. The pres- 
Sure recovery, or lack of, may also require the body shape 
surrounding the IR sensor system be reconfigured, within 


existing system constraints, to allow for ramjet operation. 
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The aerodynamic model used was that of a hemispherical body 
with an attached cylindrical body of constant diameter. The 
flow was calculated for a cylinder length equal to four 
times the nose radius. Distances involved were normalized 
to the nose radius, ra: This is the same as the radius of 
the hemisphere. 

The National Aeronautics and Space Administration Ames 
Research Center (NASA Ames) calculated the flow over a hen- 
ispherical body with an attached cylinder over the range of 
Mach numbers from 1.8 to 3.4, in steps of 0.2, at zero angle 
of attack and at Mach 3.4 with an angle of attack (a) equal 
to 10 degrees. The calculations provide flow field data 
which included shock shape, sonic line, velocity components, 
pressure, density, entropy and similar flow parameters. The 
calculations were done by computer code contained in two 
separate programs. For a further discussion of the programs 
and the various inputs see Appendix A. 

section II reports on inlet performance and the inlet 
flow field. Section III reports on ramjet performance for a 
blunt nosed centerbody for the inlet and compares the per- 


formance of a ramjet with an isentropic spike inlet. 
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II. FORMULATION OF INLET PERFORMANCE EQUATIONS 
AND INLET FLOW FIELD 


PaeewNLET DESIGN OPTIONS 

In developing a quantitative basis for comparison of var- 
10us inlet designs, several items of information are needed. 

First, information needed is the drag at the different 
Mach numbers and a correlation of drag versus Mach number. 
The blunt nosed centerbody produces a detached shock that 
stands off a distance from the body. A conical inlet pro- 
duces an attached oblique shock. 

Second, the on-design performance of the inlet design 
must be evaluated. In particular, the pressure recovery of 
the inlet is important. What happens to this recovery ratio 
when the position of the annulus is moved axially along the 
body? If the location is fixed and the radius of the inlet 
lip varied what happens to the inlet performance? A compar- 
ison must be made of the relative ranges attained by a ram- 
jet with an isentropic spike and one with a blunt nosed 
centerbody must be made. Also of concern is how sensitive 
the ramjet performance is to angle of attack variations. 

Third, off-design performance of the ramjet is of inter- 
est in several aspects. First, how does the pressure recovery 
of the inlet vary at different Mach numbers? Ramjet perfor- 


mance is sensitive to the inlet pressure recovery. Generally 
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ramjet performance tends to improve initially at higher Mach 
numbers. Second, how does this affect projectile drag and 
what effect, if any, is transmitted to the projectile range? 
Third, what is the operational region over which the inlet 
and ramjet will continue to operate? The lower performance 
limit of the inlet is that Mach number where the normal shock 
within the inlet is expelled. This condition is refered to 
as subcritical ramjet operation and can lead to a condition 
called "buzz''. The upper performance limit of the inlet is 
the maximum flight speed at which adequate pressure recovery 
can be attained. Last, how does the angle of attack affect 
off design performance? 

The desire that the ramjet powered GLM fit within the 
existing confines of the 5''/54 MK 45 gun limits the inlet 
design to that of the axisymmetric variety. In theory, 
"pop-up" inlets or scoops could be used. However, the ability 
of the inlets to survive a gun launch of 8000 g's and func- 
tion reliably increases the complexity of the round; variable 
geometry inlets are not a good use of the remaining available 
Space in the projectile. The small diameter of the projectile, 
coupled with the nose mounted laser seeker, requires intel- 
ligent. use of available space. 

The axisymmetric inlets which could be utilized are 
1) blunt nosed centerbody, 2) conical nose and variations, 

3) an unconstrained design, an isentropic spike for example, 


and 4) a cylindrical analog of the Busemann Biplane. 
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The blunt nose shown in Figure II-1l is currently in use 
on the SALGP projectile, without the annulus. As modeled by 
the simulation, the inlet consists of a hemispherical nose 
connected to a constant diameter body extending beyond the 
inlet ioe, iiemconumedis nose in Figure II-2 uses a cone of 
fixed angle to focus the shock wave at the inlet lip when 
at the design Mach number. A variation is shown in Figure 
II-3. Shown in Figure II-4 is the unconstrained design. 

The design would consist of an isentropic spike. Figure 
II-S shows the cylindrical analog of the Busemann Biplane. 
The inlet "swallows" the oblique shock wave formed by the 
nose at the design Mach number and has a body of constant 
diameter. The advantage is the elimination of external wave 


drag. 


B. CALCULATION OF STAGNATION PRESSURE RATIO OF THE INLET 

For calculation purposes, the flow around the blunt nosed 
body is divided into two segments as shown in Figure II-6. 
Region 1 is the region of flow from upstream infinity to the 
Shoulder of the centerbody. The centerbody shoulder is 
located in plane S in Figure II-6. Region 1 is characterized 
by flow in both the longitudinal and radial directions. 
Region 2 is the region of flow downstream of plane S. The 
Majority of the flow is in the longitudinal direction with a 
Delarivedy sitall radial component. The origin for coordi- 


nate systems for regions 1 and 2 is the same as illustrated 
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Figure II-1 Supersonic Inlet With Blunt Centerbody 


Figure II-2 Supersonic Inlet With Conical Spike 





Figure II-3 Supersonic Inlet (variation) 


Figure [1-4 Supersonic Inlet With Isentropic Spike 





Figure II-5 Cylindrical Analog Of Busemann Biplane Inlet 
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Figure II-6 Flow Regions Around A Blunt Nosed Body 
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in Figure 11-6. The coordinates for region 1 are X and Y, 
which are the symbols used in the NASA Ames computer printout. 
The coordinates for region 2 are R and Z, also illustrated 

in Figure II-6. Each of the regions has different reference 
values and symbols for the flow properties. Appendix A lists 
the reference values and symbols. 

In a real ramjet, the sizing of the inlet capture area, 
throat area and inlet area can dramatically vary the perfor- 
mance at different Mach numbers. The importance of the effect 
the inlet plays on performance is paramount. 

Figure II-7 shows the magnitude of pressure losses as a 
fraction of total pressure recovery. As can be clearly seen, 
the majority of loss is due to boundary layer growth. How- 
ever, as the flight Mach number increases, note that all the 
losses increase in magnitude. 

For a conical inlet, the oblique shock losses can readily 
be calculated. The overall pressure recovery value of 0.72, 
for Mach 3.0 flow, can be obtained from Figure [I-/7. 

For the blunt nosed body, BRD is tabulated in the NASA 
Ames printout. This corresponds to Pp2/Pey>? the oblique 
Shock losses, in Figure II-8. To account for the varying 
pressure recovery values within the annulus, a mass weighted 
average will be used. This concept will be discussed in 


Section II-C-3. 


28 





10 


4 
PRESSURE LOSS OUE 
O OBLIQUE SHOCKS 


ne 









° PRESSURE LOSS OUE 
bn 10 BOUNDARY LAYER 
_ GROWTH 
Us 
of 
a OVERALL 
wn INTAKE TOTAL 
us PRESSURE RECOVERY 
a. 
=. 
ts TRANSONIC SHOCK 
— AND viscous LOSSES 
SUBSONKC DIFFUSER LOSSES 
0-1 


3 b 5 5 7 
' DESIGN MACH NUMBER 
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(Reproduced from Ref. 4, Ors) 


ao 





Oblique Transonic Shock 


ee. 
a 3 


Figure II-8 Identification Of Locations For 
Description Of Pressure Recovery 
Losses In A Conical Inlet 
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iocamepresstm@ec™iarlo, Or pressure recovery, is defined 
as Ty. From Figure I[I-8, 


Pes = Pe2 «Pes 


= —=—=— = — 


(1) 
Pry Pei Pe2 


acl 
By defining 7,' as equal to P+2/P+1> the pressure recovery 


becomes 


ee ieee 
q P42 


i (2) 
The value of P,c/Py? can be determined from Figure [I-7. 
The value of Pys/Pe2 represents the accumulation of losses 
in an inlet due to the subsonic diffuser loss, Pps/Pea; 
transonic shock and viscous losses, Pra/Py3> and pressure 
loss due to boundary layer growth, P,3/Py2° The value of 


P+s/Py2 is obtained from Figure II-7 as a ratio of overall 


inlet pressure recovery to oblique shock losses. For Mach 


Su) Llow 
Pes 0.72 
Pr? Sass 0.97 — 0.7423 


Fhe value of Ta! was calculated using the flow field calcu- 
lations in the NASA Ames printout. The overall effect of the 
reduced recovery ratio is that the engine efficiency and the 


pressure ratio at the nozzle exit is changed. In turn, the 


oi 





change in the pressure ratio may cause over or under expanded 


flow from the exit nozzle [Ref. 3]. 


feeerLOW PLELD DESCRIPTION IN THE VICINITY OF A BLUNT NOSE 

For the blunt nosed centerbody at zero angle of attack, 
the flow field is axisymmetric with respect to the centerline 
of the projectile. The calculation of pressure recovery and 
ramjet performance is greatly simplified. 

1. Determination Of Capture Streamline 

Inherent to the analysis of ramjet performance is the 

determination of the mass flow rate of air at the annulus of 


the inlet, m After determining the mass flow rate, My 


L* 
the streamline coordinates will be determined at various 
points by equating the local mass flow to the mass flow at 
the lip. From these coordinates, the additive drag and 
@eenmercient of additive drag for the projectile at a parti- 
cular Mach number will be determined. Two approaches for 
calculating the streamtube points were tried and will be 
discussed here. 
a. Mass Flow Method 

The first method involves mass flow analysis. 
The mass flow through any cross section of the particular 
streamtube will equal that of the mass flow through the 
annulus. Of particular interest is the radius of the capture 


area, T,, on the bow shock wave. The air captured by this 


Buea goes into the inlet. The remainder of the air spills 
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Pa iiomualehewand atteets the performance of the projectile 
Ppa to fhe total vehicle drag with a term called addi- 
fave drag, Da: 
The mass flow at the lip is given as 


aD 
m= fr ur dA (5) 
body 


Where Py is the density, dA the elemental area and u, the 


Ik 
velocity component normal to that area. From the continuity 


equation, the mass flow at the shock front is 
° 2 
fee Ltr. uu (4) 


Refer to Figure II-9. By equating equations 3 and 4 


lip 


2 = 
ee Uw. | fr uy, 2mr, dry (cs) 
body 


By normalizing the distances involved to the nose radius and 
by dividing by po, and Via» to allow use of the NASA Ames 


printout values in region 2 , the equations become 


rou 9 ae 

ees = gee wel. 1 ar. (G3) 
Z 0 ie ae 
ice Vi co 

mm body 


Oye 
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Figure II-9 Flow Geometry Around The Blunt Nose 
Showing Bow Shock, Streamtube And 
Body Locations 
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' ' ' ; 
where u, = uy /Viy? coe ao es and r = ea The primed 
values are the tabulated values from the NASA Ames printout 
amaethe subscript refers to the applicable flow region de- 


med in Figure I1-6. Q is defined as u./V_ and equation 6 





becomes 
a 2 Pro lip ' ' ' 
> = ans fr u, T dr Cr 
xr co co 
n body 


The integrand in equation 7 was plotted and the area under 
the curve calculated for the desired lip radius. Since 


mea) l.0 inches, Q 
n 


Cc 


= 0.80174, Ores = 0.07623 at Mach 3 and 


the integral from body to lip is 0.00805, r. can be determined. 


For Mach 3 airflow, with a = Ee 2. 
: 2 
Ty) OTMBOLTTAY OT07673y «$9 00805) 


Ricmec , Ee Uo o o\. 

Table II-1 shows te fOmemach 5.0 and’a nose céenterbody radii 
equal to 1.00 inches. The maximum value of the lip radius 

is 1.854 inches. This is the size imposed by the constraints 
of the 5/54 handling system. Now the streamline locations 
aevditterent flow field points can be calculated. The compu- 
ter program used is described in Appendix B. The flow geometry 


1s shown in Figure II-9. 
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lmao et 
Capture Radius, Ts> As A Function Of Lip Radius, Ty» E Om rh 
Eauee to 1-0 Inch At Station Z=1.00 


ry (inches) ee (inches) 
0 0 OO. Z0L 
1.2000 eaoo 5 
1.5000 e705 
1.4000 . 9610 
1.5000 .1474 
1.6000 eo o 
Pee) pave ek 
1.8000 ./066 
1.8540 . 8086 
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ihemniass tlow rate across a line of constant 


angle, a, is given by the following equation 
. cs 
ee fe u, ds (8) 


where b refers to the body and cs refers to the capture stream- 
line in Figure II-9. From Figure II-10, an element of area 


on the line of constant a 1S given as 
dS = 21rdl (9) 


and dl can be related to the angle a by 


dr 


elt ern (10) 
Hence, dS in terms of dr becomes 

dS = Zur — (ie) 
From Figure II-1l1, an angle 3 is defined so that 

gB = tan t(v/u) (2, 
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Figure II-10 Magnified View Of Geometry Between 
Points 1 And 2 Of Figure II-9 
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Figure II-11 Geometry To Determine Component Of Velocity 
Normal To The Area For Integration Of Mass 
Flux 
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Hence, the normal component of the velocity becomes 


u = ue +y¥ cos(t/2-a- 8) (sy) 


Diemeche substitution of a trigonometric identity for the 


cosine term, equation 13 can be simplified to 
u = 7 u’+ V Sin (a+ 8) (14) 


The geometry of Figure [I-1l gives 


cos8 —__—_—_ (15) 


¥ ue +y 


and 


sin 8 —_____ (16) 
u“+v 


Using equations 15 and 16 and substituting into equation 14 


leaves 


5 . / ut ty2 ju sina , V_cosa } (17) 
Twa vu tv 
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tiemualuecwor atrtlow normal to the cross sectional area dS 


then becomes 
u = u sing + v cosa (18) 


Substituting equations 11 and 18 into equation 8, the follow- 


ing equation is obtained 


cs 
m = f corm) cu Sina + v cosa) /(sina) dr (19) 
b 


This becomes, after manipulation 


CS 
m = / our (ueet vv cota)dr (20) 
b 


The mass flow given by the above equation through the 
streamtube at the radial line of constant angle a must equal 
the mass flow through the nose capture area. Multiplying 

the integrand by oe and dividing both sides by the mass 


flow at the bow shock, equation 4, gives the following result 


CS 


3 ye aT u V dx 
i ~-frag’e i ge-ue (21) 
b S S n S S nh 
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In front of the projectile shoulder, equation 21 becomes 


CS 
Tr 
S ip 7 ! t 1 ' t 
a ~ * fo em cerie 7 COta) er oe) 
b 


For region 2 , which is after the shoulder, due to the dif- 
ferent reference values equation 20 becomes 
20 cs 


a 2 too ' ' ! ' ' 
CC) = On p, 1 (u, ev, Cota jidr (23) 
n 


co 60° b 





where in equations 22 and 23, the primed values are taken 
from the NASA Ames printout. In region 2 , a equals 90° and 
cotangent of a is zero. 

b. Angular Method 

The second method of streamtube calculation is 
by the angular method. A linear variation of the known flow 
field properties over small distances is assumed. The method 
Starts at the lip and progressively moves forward along the 
body determining streamtube position. 

Initially, the local velocity vector at the 
various flow field points must be plotted in polar notation. 
By averaging the velocity flow angles at two adjacent stations, 
an average angular value can be determined. As an example, 


for a lip radius of 1.2 
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Oey SCG * 84) /2 | (24) 


where 0. = local flow field velocity angle at station 3 
@, = local flow field velocity angle at station 4 


§0.,= average angle between station 3 and 4 stream 
tube locations 


Oy; the flow field angle at the lip, is determined 
from linear interpolation of the known flow field angles. By 


varying 86 with 8, known, the value of 8. is determined. 


Ao 
meas tunction of 843° 


8 = 28 = ¢ (25)) 


943 is then plotted as a function of the radius of the 
streamtube, Re in this case. 
By studying Figure II-12, the geometry of the 


field also determined Ons: 


Ry 7 Rz 
tan O43 = Ter ae J (26) 


In the above equation, all quantities with the exception of 


Rz are known. O43 is then plotted versus Rz. 
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Figure II-12 Flow Field Geometry For The Angular 
Method Of Streamtube Determination 
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The intersection of these two curves determines 
Rz and Bas» and hence 8. Gamebe determimea by interpolation. 
In this manner, the streamtube points are calculated for the 
desired lip radius. This method produces comparable results 
to the mass flow analysis but is not accurate enough pro- 
ceeding past the nore. .68 Elemprojyeetile., Additionally, the 
process is extremely slow and complicated. 

2. Determination Of Additive Drag 
The additive drag is calculated by a computer program 

discussed in Appendix C. Assuming that the points defining 
the edge of the streamtube determined by mass flow analysis 


can be connected by straight line segments, the equation for 


additive drag can be written as 
D * 2(p; - PoJA; sin@. (Zi 


where P; equals the average pressure at the data points i 

onde 1+] ; O. is the angle of the straight line segment between 

the points i and itl measured from the projectile axis; and 

A. is the surface area of the cone between the i and i+l data 

points. The geometry is shown in Figure II-13. 
inmemeoetitcient of ad@itive drag (C5442 can be de- 

termined by dividing equation 27 by q. and reference area, A 


2 i : 
r Poo (5 =e le) A, sin. 


dad YP. mM? AL 
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Figure II-13 Geometry For Additive Drag Calculation 
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The reference area, A> is taken to be the base of the five 


inch projectile. Simplification of equation (28) leads to 


C : 2s Gt) er (29) 
dad ym? ee i n 


Equation (29) is the equation programmed in Appendix C. 

The program utilizes the points generated by the pro- 
gram in Appendix B to calculate C aad’ Howenciemat 1S ampor- 
tant to scale the nose capture radius from a centerbody radius 
of one inch to the actual centerbody radius on the projectile. 
Previously r, was fixed at 1.854 inches. A ratio between 
the centerbody and lip radii is defined as x = r,/T,- 

The lip radius, Yr, , is fixed by the dimensions of the 
handling equipment, and r, can be determined as a function of 
the ratio x. Table II-2 displays the values of x, To: C aad? 
and mass flow ratio. Mass flow ratio is defined as m_/me. 

The symbols mM - and me are defined in Table II-2. 
5. Average Stagnation Pressure Ratio At The Inlet Lip 

The stagnation pressure at the annulus of the inlet 
is a function ofradius. A mass weighted average stagnation 
pressure was calculated. The average pressure recovery will 
Hewiubeiplied by a factor of P.s/P,> to obtain total pressure 


recovery as discussed in section II. 
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Table II-2 


Summary Of Radius Of Capture Streamtube, Additive Drag Coef- 
ficient, And Mass Flow Ratio 


homer = 1.654 in. 


L 

aa "7 He “n *s Caad m/m, 
tal OSL ~6855 0.6440 oO OeiZt 
ip 2 5883 -5450 0.9089 - 2956 0.240 
lo eh tO5 ~4262 1.1074 ~ 2299 Oks 57 

4 -9610 3243 eZee lv 0.471 
SS: .1474 5 5 18).0) 1.4181 s31,9 0.595 
TAG wo 5 Ss -1588 1.5473 .0928 0.697 
ey ROC OL .0906 oo! .0588 0.803 
i. 8 . 7066 .0300 © eos .0304 0.899 
1,854 . 80 86 . 0000 1.8086 20.79 0.952 


m 1s defined by equation (4); me 1s obtained from equation 
(4) with a mepidecd by ry 
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The mass weighted average of pressure recovery is 
given by 


m. 
T ee) (30) 
: Pte j=l Pto J My 


where the subscript j refers to the jth streamtube within 
ene inlet, m. is the mass flow of the jth streamtube and my. 
is the total mass flow through the entire inlet as defined 
by equation 3. This geometry is further shown in Figure 
Ti-14. If the inlet lip splits a streamtube the values for 
pressure recovery and mass flow are linearily interpolated 


to achieve the average pressure recovery. The values of 


(Pi /P too 5 were taken from the NASA Ames printout. 
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Figure [I-14 Geometry Of Streamtubes At The Inlet 
Annulus Used To Calculate Average 
Stagnation Pressure, Ta! 
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III. CALCULATION OF RAMJET PERFORMANCE 


A. PERFORMANCE EQUATIONS 

Ramjet engine performance can be quantified by the specif- 
ic thrust, specific fuel consumption, and the coefficient 
of thrust. Performance equations will be developed for each 
of the above performance criterion which include the concept 
of additive drag. 

If one dimensional flow is assumed, the thrust, F, of a 


ramjet engine is described by Netzer [Ref. 5] as 


Be: eis : » ; 
ie m.Y MVy D, A, (P, PY) (S15) 
where the subscript e refers to exit conditions and the sub- 
Scripec o refers to the inlet conditions. Term 1 is defined 
as jet or gross thrust, term 2 is the ram drag, term 3 is 
the additive drag and term 4 is the pressure thrust. If 


ideal expansion is assumed (P. = Po)» equation 1 becomes 


F=mv.-mv_- D (32) 
e O 


ee. 





the mass flow through the ramjet exit can be written as 


Mm =m, + mM, = Mm, (1 + £) (34) 
acme! 2ir patio, £, 1S generally much less than one. 
Combining equations 32 and 34 allows the total thrust to be 
written as 


F = mov, ae eV) ccrasD 3) 


a 


Paeroring Le from the first term on the right hand side gives 


F = mV (— - 1) - dD. (36) 
Combining equation 36 with the definition of Mach number and 
the speed of sound, ao» inserting an expression for combustor 
energy balance and assuming that the stagnation pressure 
through the inlet, burner and nozzle are taken as to be 
constant allows thrust to be written as a function of fuel 
air ratio, heating value of the fuel and the combustion 


tempecracture. 





F = Ve [ it CP -1 ] - D, (37) 
p to 
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The specific thrust is obtained from equation @7 by 
dividing by the mass flow rate, m. mits of m, are slug/sec 


and specific thrustis usually in units of 1b -/ (1b, /sec) . 


SF = a (38) 


The specific fuel consumption, SFC, is defined as the 
unit mass flow of fuel per hour per pound thrust. Expressed 
in the form of an equation. 


5600m¢g 


In equation (39) M has units of slug/sec. Inserting equation 





57 yields 
3600m, g 
SEC = (40) 
; h 
eee iy e+ - 1] - D 
oO =e to a 
where M = mass fuel flow rate (slugs/sec) 
g = gravitational constant 
h = heating value of fuel (BTU/1b) 
f = fuel air ratio 
. = specific heat capacity (BTU/1b °R) 
iy = inlet stagnation temperature 


So 





mie tnuist: COerricient is defined as 





_ i 
C = A (41) 


where do 1s the dynamic pressure and A, is the reference 
area. the reference area is the base area of the five inch 
projectile. For a ramjet, mass flow at the inlet is given 


as 
m = AAV p (42) 


Combining equations (37), (41) and (42) gives 


2A, ac dD. 
o = —— [/v¥1l1 + - lj - (43) 
rE CO to Toh 








The preceding equations were derived for a ramjet engine 
without internal losses. Losses in a ramjet engine that 
Should be accounted for are 


a) loss of stagnation pressure in the diffuser, T 3 
b) loss of stagnation pressure in the burner, TH 
c) loss of stagnation pressure in the nozzle, ie 
d) combustion efficiency, Ny 
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These terms result in the generation of a common term, 
Si which does affect the calculation of ramjet performance. 


5, is defined as 





l y-l W/Z 
( \igaee = 1 
, THT aT nth 
aan cl ome hh ht ct eg 
ae pec° 


The square root appearing inequations 68), @0) and @3 is replaced 
by S. to give ramjet performance with internal losses and with 


an additive drag term. The performance equations becomes 





SF = [mM a, (S$, - 1) - DJ]/mg (45) 
3600m -g 
2A, Do 
Ce i oa (S.. 7 1) = q A (47) 
ue OT 


The computer program developed by Fuhs [Ref. 6] was 
modified to compute the additive drag term. Holding the values 
of T > TT, and n, constant, the performance of a ramjet engine 
with a blunt centerbody was predicted. The performance of 
a ramjet with a conical spike inlet as predicted by the 


program serves as a baseline measurement of performance. 
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Teo oelSsetON OF RESULTS 


A. PERFORMANCE CRITERION 
Performance of a ramjet engine can be quantified in 
different ways. For purposes of this report, the performance 
was compared to a ramjet baseline configuration uSing a coni- 
cal inlet spike; the inlet capture area is 0.0123 £t? [Ref. 7]. 
The performance of a ramjet with the blunt centerbody 
configuration was calculated for a range of fuel/air ratios. 
Performance is specified on the basis of combusiton exit 
temperature (T 3) specific fuel consumption and excess 
thrust. Of special interest is performance when Cr equals 
c 


The value of Cy» which was obtained from White [Ref. 8], 


D* 
was 0.349. 

The combustor exit temperature is calculated by the ram- 
jet performance program as an output variable. Current tech- 
nology limits the steady combustor exit temperature to about 
4400°R. The value chosen for analysis purposes is 4422°R. 
Above this temperature, the combustor and exit nozzle will 
melt, if run for a continuous period of time. The value of 
T.3 is generally less than T.3 at the stoichiometric fuel/air 
ratio for the hydrocarbon. 


roreeive urinal seconds of flight, it 1s possible to boost 


the fuel/air ratio to the stoichiometric value. The boost 
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in fuel/air ratio will increase the amount of excess thrust 
available for terminal maneuvering. Since the GLM will 
explode at flight termination, the combustor exit temperature 
memot little importance in terminal flight. 

A further basis for performance comparison is the effect 
on relative detection range of seeker aperture. Seeker 
aperture is equal to mr”. 

Specific fuel consumption (SFC) was chosen as a compari- 
son standard because it provides a good measure of ramjet fuel 
efficiency. Excess thrust is used as a measure of the ability 
of the guided projectile to maneuver in the terminal phase of 
the encounter. Regardless of the seeker accuracy, lack of 
thrust above the amount of flight drag will result in the 


inability to pursue a maneuvering target. Excess thrust 


coefficient, Cees is calculated as 


Cee Tez) = Sg (TE) - Sp (48) 


where Cy 1s airframe drag and C -(Te3) is the thrust coeffi- 
cient for Th3: Two values of T,3 acemOMmErMcereste. One 15 
the maximum allowed combustor exit temperature for steady 
operation which is 4422°R. The other is T.3 for stoichiomet- 


ric combustion which is 5968°R; stoichiometric fuel/air 


ratio is 0.06. 


>/ 
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B. RESULTS OF CALCULATIONS 

Figure IV-1 plots T,, versus ry/t, ey oot Cy = Ce. 
As can be clearly seen, as the ratio of ye increases, 
or as the centerbody shrinks, T.3 decreases. T Ls reaches 
the maximum sustained operating value of 4422°R when ry, /T, 
omer OS ; r,/T, equal to 1.63 is the minimum value for sus- 
tained operation. For stoichiometric combustion (f=0.06), 
t/t, could possibly be reduced to a value as small as 1.50. 

From the minimum value of r,/T, obtained in Figure IV-l, 
Preure I1V-Z is entered. Figure IV-2 is a plot of excess 
thrust coefficient versus ey and of stoichiometric excess 
Gheust coetficient versus r,/T- Stoichiometric excess 
thrust is defined as Ce at stoichiometric conditions (f=0.060) 
minus Cy: 

Using a value of Eee equal to 1.63, the ramjet will 
have excess maneuvering thrust only at values of rj /T, greater 
than 1.63. Use of stoichiometric fuel/air mixture, during 
the last seconds of the encounter with a maneuvering target, 
results in excess thrust being generated for r,/T, greater 
enamel. S51. 

Figure IV-3 is a graph of the specific fuel comsumption 
(SFC) versus r,/T,- SFC decreases as ee increases. The 


baseline ramjet has a SFC equal to 2.02 for C.,, equal to Ce. 


D 
For ae equal to 1.63, the ramjet with a blunt nose has 


65% higher SFC. 
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Figure I[VY-1 Combustor Exit Temperature As a Function 


Of ae At Ce = Cy = 0.349 
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Specific Fuel Consumption As A Function 
Of r,/T, For Ramjet With Blunt Centerbody 
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Figure IV-4 is a plot of relative detection range of an 

IR seeker for various r,/T,- The signal-to-noise ratio of 
the IR seeker and, therefore, the detection range is propor- 
tional to the area of the lense, The lense area is directly 
related to the square of the lense radius. As Ieee in- 
creases, the detection range decreases as (r,/t,)°. 

Mass flow ratio was defined in the text near Table II-3. 
The ratio is interpreted as the fraction of the mass flow 
actually captured, m., relative to mass flow into an area 
- 
of r,/T,: Also Figure IV-5 has a plot of capture streamtube 


Figure IV-5 has mass flow ratio plotted as a function 


fadeus ratio, ro/Ty > a5) a GunGeion of r,/T,- For the mini- 
mum value of r,/T, equal to 1.63, the capture radius ratio 
is 0.84. Likewise, the mass flow ratio is 0.70. 

Figure IV-6 shows the variation of additive drag coeffi- 
cient with r,/T,- At the minimum value of pe Ee Of 1 2G, 
Coad Powweco. The airframe drag coefficient is 0.349. Hence, 


C ms 24% of Cy: To understand the behavior of Caad as a 


dad 
fumetion of r,/T, refer to Figure IV-7. For r,/T, equal to 
1.1, which is illustrated in Figure IV-7(a), the capture 
streamline between bow shock and lip of the annulus is steep. 
Also, the pressure is large near the axis. Consequently, 
Cea is large for small r,/T,- 

The greater the ratio of r,/T, becomes the better ramjet 


performance becomes. However, to achieve even marginal 


performance at r,/T, = 1.63, anwinrealistic penalty im the 
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Figure IV-4 Relative Detection Range As A Function Of eae 
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form of a 63% loss of relative detection range must be paid. 
Lower values of se) Se are unable to be obtained in a thrust 
equals drag configuration due to high values of additive 
drag and poor Ta 

Table IV-1 contrasts ramjets with the spike inlet and 
two blunt nose inlets. All three ramjets have identical 
thrust coefficients equal to airframe drag coefficient. Due 
to large specific thrust, the ramjet with spike inlet re- 
quires considerably less mass flow rate. The ramjet with 


spike inlet has superior specific fuel consumption. Excess 


thrust is comparable for the three ramjets. 
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Table VI-1 
Comparison Of Spike Inlet Versus Two Blunt Nose Inlets 


eee 
Spike 
Inlet Blunt Nose Inlet 





Parameter 


Airframe drag 


coefficient, Cy 0.349 0.349 0.349 
Additive drag 

coefficient, C 0 0.018 0.085 

dad 

Thrust coefficient, Ce 0.349 0.349 0.349 
Combustor exit temperature 

a T.3 AeA 3284 4422 
Mass flow rate, 1b, /sec 9.5 1355 13.410 
poeeriic Thrust, 

1b -/ (1b, /sec) G7. 2 Ses 46.8 


specific fuel consumption, 
(1b /hr)/1b, Dee Zos SS 


Maximum thrust coefficient 
ee, at stoichiometric 
fe air ratio 0.505 0.798 0.509 


Combustor exit tempera- 


ture, T.3> °R 5698 5698 5698 
Relative detection 
range Oro 3 OR29 0.37 


* Best ramjet is obtained for r,/r_ equal to 1.854 at least 
for the range of r,/r investigated here. Best IR 
detection capability Sccurs for eae equal 1.63. 
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V. CONCLUSIONS 


Conversion of the Navy SALGP from a rocket to a ramjet 
severely degrades detection capability, i an axisymmetric 
inlet is used. Note the values of 0.33, 0.29, and 0.37 for 
relative detection range from Table IV-l1. 

The ramjets with blunt nose inlets suffer in both speci- 
fic thrust and specific fuel consumption. The poor perfor- 
mance is due to two causes. First, the high value of oa 
and second, poor pressure recovery of the inlet leading to 
poor performance. 

Surprisingly, the excess thrust coefficient of ramjets 
with blunt noses is very competitive with the ramjet with 
Spike inlet. 

In view of the foregoing conclusions, pop-out or retrac- 
table scoop inlets appear much more attractive in spite of 
structural constraints and space considerations. Use of an 
axisymmetric inlet located in the nose of a 5-inch projectile 


1S a severe design compromise. 
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APPENDIX A 
NASA AMES COMPUTER PROGRAM DESCRIPTION 


I. INTRODUCTION 

Two different computer programs were used to model the 
flow around the blunt nosed body. The coordinate system 
used in the thesis and the NASA Ames computer program is 
shown in Figure [I-6. 

The first, IMPLCBO, is a modified version of AXI-BLUNT; 
see Kutler, Chakravarthy, and Lombard [Ref. 9]. Programmed 
to predict the supersonic flow over a three dimensional 
body, the algorithm predicts shock shape and location as well 
as flow parameters at equally spaced interior points between 
the body and the shock. 

The second program, OGIVE, solves for the supersonic flow 
parameters around a three dimensional wing body configuration; 
refer to Kutler, Reinhardt, and Warming [Ref. 10]. Originally 
designed for flow prediction around a delta winged spacecraft, 
the program will provide flow field parameters at various radii 
from the centerbody as well as at various locations along 


the body. 


II. PROGRAM DESCRIPTION 
The IMPLCBO program uses an unsteady, implicit numerical 
procedure to determine the supersonic flow around the body. 


Further program description is contained in Reference 9. 
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The program output is referenced to various free stream 
values. These values are shown in Table A-l. 

The program utilizes the body configuration and predicted 
shock points as boundary values and solves the inviscid 
Euler or Navier-Stokes equations for the interior points, 
after a coordinate transform. 

The different angular relationships at the shock front 
are shown in Figure A-1l. 

The OGIVE program utilizes the flow field parameters 
determined by IMPLCBO, at 9 = 90° and X = 1.00, as an initial 
starting point and determines flow field values along the 
cylinder body by a finite difference algorithm. The program 
is further described in Reference 10. The program reference 
values are tabulated in Table A-1l. 

The program describes the body shapes by analytical approx- 
imations. These approximations consist of coordinate posi- 
peomemandecne Slope of the body contour as a function of the 


distance along the body axis. 
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Table A-l 


Reference Values For OGIVE And IMPLCBO Programs 


IMPLCBO (region 1) 
output variable 
P/PINF 
RHO/RINF 


U/QINF 
V/QINF 


OGIVE (region 2) 


R 
fe 


RHO 


QINF 


reference 
Po 

Po | 
q, = vyYM 
q, = VYM 


Taek. 0 
nh 

Pro > 
Gai), ve 9 
qn 
a ; 
Gn 
Pn ; 


fz 


es Gor 


value 


Poof Pew 0 -0272237 at 
Mine iO 

Pig/ P47 0.0762263 at 
Peron 

2yl/2 

Z component of flow 

VielOc il en, 

R component of flow 

Weal OG 1 ty 

Q component of flow 

velocity 

length along 

cylinder body 





Bow Shock 


B 


: 
\ 


Body 


Figure A-1 Angular Relationship At The 
SHOCK Front For 6, 8, 06 
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APPENDIX B 
MASS CONTINUITY METHOD OF DETERMINING STREAMLINE SHAPE 


A computer program has been written for the HP 9830 to 
Calculate the capture streamline. The program is designed 
to reduce the output data obtained from the NASA Ames com- 
puter output and compute the streamtube position at various 
points in the flow field around the blunt nosed body. 

The streamtube position is determined by comparing the 
mass flow of air entering the bow shock to the integral of 
mass flux. The streamtube radius at the bow shock, Tos 1s 
determined by comparing the mass flow at the inlet to mass 
flow at the nose. 

Miewvalue of the integral in equation (22) or (23) is 
determined by assuming a linear change in pressure, density 
and velocity between successive data points. The program 
interpolates to find the distance from the axis that the 
streamtube must be to satisfy the mass flow requirements of 
the inlet. Values of X and Y as shown in Figure II-6 are 
calculated, and an interpolated pressure is determined. 
The value of the integral in equation (22) and (23) is cal- 
culated using the rectangular method. Equations (22) and 


(23) are developed in Chapter III of this thesis. 
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Table B-l 


Variables For Capture Streamtube Computer Program 


Symbol for 
HP 9830 


AO 


A7 
C 


Cl 
GZ 
IP 


G 
I 
P7 
S7 
W 
Q9 
RO 
S 


X7 


Ne 


Definition 


value of integral in equation (22) oh ae AD 
between point i and itl 


remaining area between tabulated values 


squared radius of streamtube at shock front 
(square inches) 


constant for integral when in region I 
constant for integral when in region II 


percentage of distance Y7 is between Y(I) 
and Y(I+1) 


dummy variable for area calculation 
indexing variable 

imceerpolated value of pressure X7,Y7 
slope of area curve 

dummy variable used in calculation of C 
value of Q., equal to (1+ (Gin 12 
value of p,/ Ome at M, 


Summation of integral area to the ith point 
(inches squared) 


interpolated value of X, determined from Y7 
(inches) 


interpolated value of Y (inches) 


Input variables 


A 


ee 
z 


angular value of data points above projectile 
centerline (degrees) 


number of data points in data file 


location of data points in relation to the 
projectile shoulder. 1= in front of shoulder 
2= on the body cylinder 
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EEE 


Table B-l con't 


Symbol for. Definition 
fe 9830 


Dimensioned Variables 


Gi} iieecral Value between 1 and itl] 

Iie } integral value at ith point 

alee ) pressure ratio at ith point 

Ret ) Geiisuty fatio a2t rth point 

UCI) component of flow velocity along the body 


in the Z-direction 


v1) component of flow velocity along the body 
Pieetenr direceron 


at) distance from the body centerline to the 
data point 
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APPENDIX C 


ADDITIVE DRAG COMPUTATION 


The program is designed to calculate the additive drag 
and coefficient of additive drag for an inlet. The reference 
area is the projectile base area. The program is written in 
Perc tor the HP 9830 desk top calculator. 

The input data consists of the point position (X and Y) 
and the static pressure at the point. The three variables 
X, Y, and p are determined by first running the program 
described in Appendix B of this thesis. 

Equations (27) through (30) used in this calculation are 
developed in Section III of this thesis. 

The computer symbols are listed and defined in Table C-l 


and the program listing is in Table C-2. 
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Table C-l 


feeraples For Coefficient Of Additive Drag Computer Program 


Variables 


A 
C 


ce 


Pi 


° 
ol 


W 


Definition 


Reference area (square inches) 


Coefficient of additive drag between I and I+l 
data point (dimensionless) 


total coefficient of additive drag (dimensionless) 
total drag of body 


distance along body, Z, between data points 
(inches) 


counting variable 

counting variable 

Ratio of circle circumference to diameter 
ober at inpeie M 

Distance from centerline, Y (inches) 

Slant distance between data points (inches) 


Surface area of cone between data points 
(square inches) 


Incremental drag between data points (1b) 


Input variables 


M 
in 


Mach number (dimensionless) 
Number of data points 


Dimensioned variables 


Patel ) 
(a ) 
Rel ) 


X(T) 
(1) 


Angle of streamtube between data points (degrees) 
Pressure at Ith data point (dimensionless) 


Average pressure between data points 
(dimensionless) 


X position of Ith data point (inches) 
Y position of Ith data point (inches) 
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